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ABSTRACT 



Aims. This work aims to disentangle the morphological, kinematic, and chemical components of the nebula Ml -67 to shed light on 
its process of formation around the central Wolf-Rayet (WR) star WR124. 

Methods. We have carried out integral field spectroscopy observations over two regions of Ml -67, covering most of the nebula 
in the optical range. Maps of electron density, line ratios, and radial velocity were created to perform a detailed analysis of the two- 
dimensional structure. We studied the physical and chemical properties by means of integrated spectra selected over the whole nebula. 
Photoionization models were performed to confirm the empirical chemical results theoretically. In addition, we obtained and analysed 
infrared spectroscopic data and the MIPS 24//m image of Ml -67 from Spitzer. 

Results. We find that the ionized gas of Ml -67 is condensed in knots aligned in a preferred axis along the NE-SW direction, like a 
bipolar structure. Both electron density and radial velocity decrease in this direction when moving away from the central star. From 
the derived electron temperature, T e ~8200 K, we have estimated chemical abundances, obtaining that nitrogen appears strongly 
enriched and oxygen depleted. From the last two results, we infer that this bipolarity is the consequence of an ejection of an evolved 
stage of WR124 with material processed in the CNO cycle. Furthermore, we find two regions placed outside of the bipolar structure 
with different spectral and chemical properties. The infrared study has revealed that the bipolar axis is composed of ionized gas with 
a low ionization degree that is well mixed with warm dust and of a spherical bubble surrounding the ejection at 24//m. 
Taking the evolution of a 60 Mo star and the temporal scale of the bipolar ejection into account , we propose that the observed gas 
was ejected during an eruption in the luminous blue variable. The star has entered the WR phase recently without apparent signs of 
interaction between WR- winds and interstellar material. 

Key words. ISM: bubbles - ISM: abundances - ISM: kinematics and dynamics - ISM: individual: Ml-67 - Stars: Wolf-Rayet 



1. Introduction 

WR124 (=BAC209) is a Galactic massive star characterize d by 
a very high heliocentric recession velocity of -175 km s" 1 (Solf 
|& Carsenty||1982| ), and it is regard ed to be among the fastest 
moving mas sive stars in the Galaxy ( Moffat et al.|1982| ). It was 
classified by Merrill (1938) as a nitrogen-sequence Wolf-Rayet 
star (WN) and later as Population I WN8 star ( |Smith & Afler| 
[1969). 



1965). The study of nebulae around WRs gives us clues to the 
mass-loss history of massive stars, as well as to the chemical 
enrichment of the interstellar medium (ISM). 



WR stars are thought to be a late stage in the evolution of 
stars more massive than 25 Mq and they are characterized by 
significant stellar winds with high mass-loss rate and terminal 
velocity. Many WRs are surrounded by nebular emission, some 
of which are members of a class of object called ring nebulae. 
The structure of this type of nebula is attributed to a continual 
process of mass loss from the exciting WR star, which sweeps 
the surrounding interestellar gas into a shell (Johnson & Hogg 



* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck-Institut fur Astronomie and the Instituto de Astrofisica de 
Andalucia (CSIC). 
** e-mail: alba@iaa.es 



Ml-67 (=Sh2-80, Sharpless| ^959] ) is a bright nebula sur- 
rounding WR124, and it shows a clumpy and irregular distri- 
bution of gas that is mostly condensed in bright knots and fila- 
ments. It was first reported by Minkowski (1946) during an Ha 
objective prism survey. Classification of the nebula and its dis- 
tance have been subjects of debate in the past years. Although it 
was first considered an Hn region (Sharple s~s|1959| ), the classifi- 
cation of Ml-67 has been alternating between a planetary neb- 
ula (PN) and a ring nebula. Bertola (1964) adopted a distance 
of 0.9 kpc and suggested that Ml-67 might be a PN since both 
star and nebula have high radial velocity. Studies from optical, 
infrared, and radio data by |Cohen & Barlow| ( |1975| ) prompted 
its classification as ring neb ula around a WR, with a distanc e of 
4.33 kpc (in agreement with Pismis & Recillas-Cruz| ( |1979] ) es- 
timations). Nevertheless, |van der Hucht et al.| ( [1985| ) supported 
the PN status based on the energy distribution in th e far infrared. 
The issue was settled by Esteb an et al.| (T991) and C rawford &| 
Barlow (1991). The detailed abundance analysis of the nebula by 
Esteban et al.| ( |1991| ) revealed nitrogen enhancement and oxy- 
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Table 1. Main physical parameters of WR 124 and Ml -67. 



Object Parameter 



WR124 



Ml-67 



Value 



Reference 



(a,6) (J2000) 


(19:11:30.88, +16:51:38.16) 


Spectral type 


WN8 


Voo (km s _1 ) 


710 


Tcs (kK) 


44.7 


Distance (kpc ) 


4-5 


R G (kpc ) 


8-10 


v he i (kms -1 ) 


175 


M v (mag) 


-7.22 


E h _ v (mag) 


1.08 


Ha diameter (arc sec) 


110-120 


v he i (kms -1 ) 


150-185 


v exp (kms -1 ) 


46 


M ionized (Mq) 


1.73 



Perryman & ES A 



Smith & Aller 




van der Hucht 
Hamann et al 
Crawford & Barlow (1991) 
Esteban et al. (1992) 
Solf&Carsenty(1982) 
Hamann et al. (2006) 
Hamann et al. (2006) 



Grosdidier et al.]<|1998| > 



Chu & Treffers (1981) 



Sirianni et aHl9 98) 
|GrosdidieretaL](1998 ) 



gen deficiency, which is typical of material ejected in a previous 
evolutionary phase, and pointed to a progenitor more massive 



than those usually associated with PN central stars. Crawford & 



Barlow (1991) estimated a distance between 4 kpc and 5 kpc us- 



ing the interstellar NaiD 2 absorption spectrum of the star, ruling 
out the PN nature. Recently, Marchen ko et al.| ( [2010| ) have used 
a comprehensive model of the nebular expansion to estimate a 
distance of 3.35 kpc. Currently, Ml-67 is classified as an ejected 
type WR ring-nebula. 

Although Ml-67 shows an apparent spherical symmetry, 
ground-bas ed coronographic images revealed a bipolar structure 
( Nota et al. |1995a| ). The emission lines seems to be caused by 
condensations of gas in clumps and radial filaments ( [Grosdidier 



|etal.|1998] ). One of the most striking characteristics of the nebula 
is the virtual absence of optical oxygen emission lines (Barker 
[1978; Esteba ri^raT|T99TT ). Nevertheless, [Chu & T reffers (198T) 
reported a bright spot of [Om]/i5007A 15" to the NE of the cen- 
tral star. Spectroscopic investigations of the physical conditions 
and abundances of the nebular shell have shown that the ionized 
gas is nitrogen-enriched and oxygen-depleted, suggesting that O 
has been processed into N mainly via the ON cycle (Esteban et 
aT7||1991) ). This implies that Ml-67 is almost completely com- 



posed of stellar material that is poorly mixed with the surround- 
ing ISM. The long-slit spectroscopy of Ml-67 established that 



the bulk of the nebula is expanding at y exp =42-46 km s 1 (|Solf & 



Carsenty|1982[ |Sirianni et al.|1998) and |Chu & Tr effers 



(Soli& 



confirmed the high heliocentric velocity of the nebula Vh e i=150- 
185 km s" 1 , which is comparable to the velocity of the star. The 
main parameters of the central star, WR124, and the nebula, Ml- 
67, are summarized in Table [T] 

Many studies have tried to disentangle the geome try and 
dynamics of Ml-67 and its interaction with WR124. |Solf & 
|Carsenty| ( |1982| ) proposed a simple expanding "empty" shell 
with condensation of stellar material that was ejected by the 
high- velocity parent star; indeed, the leading edge of the shell is 
considerably brighter than the trailing part. Siri anni et al.| ( |1998] ) 
found two components in the environment of the central star and 
interpreted them as the consequence of two different events in 
the past: a spherical hollow shell of 92" in diameter expanding 
at 46 km s" 1 and a bipolar outflow with a semi-dimension of 48" 
and a velocity of 88 km s" 1 with respect to the expansion centre. 
On the other hand, some authors explained the asymmetr y as 
the result of a possible low-mass companion for WR124 ( [Chu 



|& Treffers|198T]|Moffat et al.|1982| ). 

Despite the important findings of the last years, some 
relevant aspects of the evolution and formation of the ring 
nebula associated with WR124 remain unknown. In particular, 
a 2D study of the ionization structure of the nebula covering 
all the morphologies and/or the structural components can shed 
light on the formation process of the nebula from the ejecta 
of the central star. The late spectral type of the WR ionizing 
star (WN8) is also very remarkable to study, as well as the 
degree of homogeneity in the chemical composition of its ejecta. 

To do this, we included Ml-67 in our programme of integral 
field spectroscopy (IFS) observations to compare the 2D struc- 
ture with the integrated properties of certain selected areas and 
with models of WR evolution. The paper is organized as fol- 
lows. First, we describe the observations and data reduction in 
Sect. [2] Then, we present the 2D results for morphology, ioniza- 
tion structure, and kinematic in Sect. [3] In Sect. [4] we show the 
physical conditions and chemical abundances of eight selected 
areas. We perform a study of Ml-67 in the mid-infrared range 
by analysing the IRS spectrum and the 24/mi MIPS image from 
Spitzer in Sect. [5] In Sect. [6] we discuss the chemical composi- 
tion of Ml-67, the observed structure, and its relation with the 
evolution of the central WR star. Finally, a summary of the main 
conclusion is given in Sect. [7] 

2. Observation and data reduction 

The observations were carried out on July 5, 2005 using the 
Potsdam Multi- Aperture Spectrograph instrument (PMAS) 
([Roth et al 112005] ) in PPAK mode (PMAS fibre Package, [Kelz 
et al. 2006) at the 3.5 m telescope of the Centro Astronomico 
Hispano Aleman (CAHA) at the observatory of Calar Alto 
(Almeria, Spain). The PPAK fibre bundle consists of 382 
fibres with a diameter of 2.7 arcsec. The 331 science fibres are 
concentrated in an hexagonal bundle covering a field of view 
(FoV) of 74" x65". The surrounding sky is sampled by 36 fibres 
distributed in six bundles located following a circle at about 
90" from the cent re. There are 15 fi bres for calibration purposes 
too (see Fig. 5 in jKelz et aL]|2006| ). We used the V300 grating, 
covering the spectral range from 3660 to 7040 A with a disper- 
sion of 1.67 A/pix, giving a spectral resolution of FWHM-8.7 
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Fig. 1. Narrow-band image of Ml -67 in Ha+continuum taken with the Wide Field Camera at the Isaac Newton Telescope. North is 
up and east left. Red hexagons show the two zones of our IFS observations: Edge to the NE (left) and Centre to the SW (right). 



A (R = A/AA ~ 660) at 5577A. The weather was photometric 
throughout the observations with the typical seeing subarsecond. 

To choose the regions of Ml -67 to be mapped, we resorted 
to the narrow-band images observed by our group at the Isaac 
Newton Telescope (INT) with the Wide Field Camera (WFC). 
The first PPAK pointing (called Centre) was centred on the WR 
star and it covers almost the whole nebula. The second zone 
(called Edge) was selected to study the NE edge of the object 
containing nebular emission and surrounding medium. Both re- 
gions can be seen in Fig.[T] Table [2] shows the observational log 
for Ml -67. 

Bias frames, continuum, arcs, and one spectrophotometric 
standard star (Hz 44) were also acquired during the observations. 



Th e data were reduced using the R3D software (Sanchez 
2006| in comb ination with IRAffJ and the Euro3D packages 



(Sanchez 2004). The reduction consisted of the standard steps 
for fibre-fed IFS observations. 

At first, a master bias frame was created and subtracted from 
all the images. The different exposures taken at the same po- 
sition on the sky were combined to reject cosmic rays using 
IRAF routines. A trace mask was generated from a continuum- 
illuminated exposure, identifying the location of each spectrum 
on the detector along the dispersion axis. Then, each spectrum 
was extracted from the science and standard star frames, co- 
adding the flux within an aperture of five pixels at the location 
of the spectral peak in the raw data using the tracing informa- 
tion, and storing it in a 2D image called row-stacked spectrum 

1 The Image Reduction and Analysis Facility IRAF is distributed by 
the National Optical Astronomy Observatories, which are operated by 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. Website: 
http ://iraf . noao . edu/ 



(RSS) (Sanchez 2004). We checked that the contamination from 
flux coming from adjacent fibres using this aperture was negli- 
gible (Bershady et al. 2004; Sanchez 2006). For a given aperture 
and FWHM ~ (0.5 x aperture), we found a level of cross-talk 
that was always <10%. This seems to be an acceptable compro- 
mise between maximizing the recovered flux and minimizing the 
cross-talk. 

Distortion and dispersion solutions were obtained using a 
He calibration-lamp exposure and applied to the science data 
to perform the wavelength calibration. The accuracy achieved 
was better than ~0.1 A (rms) for the arc exposures. Corrections 
to minimize the differences between fibre-to-fibre transmission 
throughput were also applied, creating a fibre flat from the ex- 
posure of a continuum source. Observations of the spectropho- 
tometric standard star Hz 44 were used to perform the flux cali- 
bration. 

The sky emission was determined using the science data 
spectra obtained throughout the additional fibres for sampling 
the sky. As explained above, the second pointing was made at the 
edge of Ml -67, and some of its sky bundles are located within 
an area containing signals from the nebula. We inspected the 
36 sky-fibres of each pointing, selecting those that did not show 
nebular emission. The spectra of all the selected fibres were com- 
bined with a mean in a single spectrum, and a 2D spectrum was 
created by copying the combined spectrum in each fibre. These 
sky spectra were then subtracted from every science spectrum, 
each pointing with its own sky. 

Finally, considering the wavelength range and the airmass of 
the observations, and using Filippenko ( 1982), we estimated that 
offsets due to the differential atmospheric refraction (DAR) were 
always smaller than one third of the fibre diameter. Correction 
for DAR was not necessary in our data. 
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Table 2. Ml-67 PPAK observational log. 



Zone 


Coordinates (J2000) 
(ar,<5) 


Grating 


Spectral range 

(A) 


Exp. time 

(s) 


Airmass 


Date 


Centre 
Edge 


(19:11:30.9, +16:51:39.2) 
(19:12:14.8, +16:52:12.9) 


V300 
V300 


3640-7040 
3640-7040 


3x30 
3x450 


1.08 
1.07 


July, 5, 2005 
July, 5, 2005 



3. Two-dimensional analysis 

To perform a detailed analysis of the 2D structure of the nebula, 
we built easy-to-use datacubes of the two reduced pointings 
with three dimensions: two spatials and one spectral. We studied 
all the interpolation routines included in the E3D software and 
verified which conserves the flux and the apparent morphology 
observed in the spaxels. Finally, we generated our cubes with the 
linear Delaunay interpolation and a pixel size of 1 .5x1 .5 arcsec 2 . 

In a first attempt to understand the morphology of the two 
observed zones, we extracted images by sliding the cubes at dif- 
ferent wavelength ranges. They are presented in Fig. [2] on loga- 
rithmic scale. 

In the 5006A-5014A range, which includes the [Om]/i5007A 
line, no significant extended emission can be observed in both re- 
gions, supporting previous studies that revealed no oxygen emis- 
sion in Ml-67. Several spots appear in the FoV (including the 
central WR124 star) with fluxes lower than ~ 10" 17 erg cm -2 s" 1 , 
but we checked that their nature was not nebular. They probably 
are stars in o ur line of sight. W e gave special attention to the spot 
described by |Chu & Treffers| ( ^98T] ) at 15" to the NE of the star. 
Although we can observe some emission, we cannot confirm that 
it comes from the nebula. A more detailed analysis of this spot 
is performed in Sect.[4]by means of the integrated spectrum. 

As for other lines from the central pointing maps (Ha, [Nn], 
and [Sn]), most of the emission seems to be concentrated in 
at least five knots distributed in the NE-SW direction without 
counterpart in the [Om]/i5007A image. In addition, two regions 
with very faint surface brightness (or even no emission) can 
be seen at the opposite sides (NW and SE). This orientation 
agr ees with the bipolar str ucture observed by Nota et al. (1995a) 
and Sirianni et al. ( 1998| ) with coronographic studies. Ha and 
[Sn] emission shows a discontinuity in the edge pointing of the 
nebula with higher surface brightness in its SW area. When 
we move to the NE, the emission decreases until it disappears. 
The purple coloured area reaches non-negligible emission up to 
Ha ~ 10" 16 erg cm -2 s" 1 per pixel (lpixel=2.25 arcsec 2 ). 



3.1. 2D study of the emission-line maps 

Maps were created from cubes by fitting the emission lines 
in each spatial element fol lowing the methodology presented 
~( [20T2| ). 



CENTRE 



in 
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Basically, we performed a 
Gaussian fit to the emission lines using our own routine, which 
returns maps of the flux, centre, and FWHM among other prop- 
erties. To prevent contamination by low signal-to-noise (S/N) 
data, we masked out all pixels with S/N lower than 5. During 
the creation of the S/N masks, we visually inspected all the pix- 
els, rejecting interactively those with non-nebular spectra or with 
contamination from the central WR and other field stars. 

For both regions (centre and edge pointings), maps of 
parameters from the Gaussian fitting were generated for 




Fig. 2. Interpolated images of Ml-67 of the two observed re- 
gions: In the left column the edge pointing and in the right the 
central one. In each row we represent the flux (including con- 
tinuum) integrated in a wavelength range. Top: range 5006A- 
5014A including [OinU5007A. Middle: range 6562A-6590A in- 
cluding Ha and [Nn]/16584A. Bottom: range 6729A-6737A in- 
cluding [Sn]/L167 17,6731 A. All the maps are represented on log- 
arithmic scales with units of log(erg cm -2 s" 1 ). The size of the 



hexagon side is 38 r 
left (see Fig. [I]). 



In all the maps, north is up and east to the 



seven emission lines: Hy, H/3, Ha, [Nn] /U6548,65 84 A, and 
[Sii]/L167 17,6731 A. Although our spectral range includes the 
[Oii]/Li3726,3728A lines, their automatic fit was not considered 
because these lines are faint and placed at the edge of the CCD, 
where the distortion correction bended and deformed them. 

All the emission line maps were reddening-corrected using 
the reddening coefficient, c(H/3), map (each pointing was cor- 
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Fig. 3. Spatial structure of the derived c(UJ3) maps and their cor- 
responding statistical frequency distributions with a binning of 
0.1. On the left the edge pointing and on the right the central 
one. Orientation and sizes of maps are as in Fig. [2] 



rected with its own c(H/3) map). To determine this coefficient 

we resorted to the Ha/H/3 line ratio. We analysed the maps of 

the three B aimer lines detected in our wavelength spectral range 

(Hy, HJ3, and Ha) and decided to discard the Hy flux because 

the S/N was lower than 5 in -20% of the spaxels. However, 

we checked that both derivations were consistent in the spaxels 

with good S/N. We used an intrinsic Balmer emission line ratio ( |1975| ) agree with our values. 

of Ha/Hj3 =3.03 obtained from the public software of Stor ey& 



Fig. 4. Electron density, n e , maps derived from the 
[Sn]/U67 17/6731 line ratios in units of cm -3 . Orientations 
and sizes as in Fig. 2] On the bottom the statistical frequency 
distributions with a binning of 100 cm -3 . The edge pointing is 
on the left, and the central on the right. The white lines across 
the maps (from NE to SW) represent the direction along which 
the cuts were extracted to study the radial variation of n e . See 
text for details. 



aLl(pMi) (-1. 35), whereas the A, 
Tref fers| ( |198l) and E(B 



V) 



j - 3.8 o btained by |Chu & 
1.35 from jCohen & Ba rlow 



Hummer ( 1995 ), assuming Cas e B recombination with an elec- 
tron density of n t - 1000 cm -3 (Esteban et al. 1991 ) and an elec- 



tron temperat ure of T Q - 70 00 K (the mean value between the 
estimations of |Barker|[T978l (~7500 K) and Esteb an et aL][T99Tl 
(-6000 K)). Statistical frequency distributions of the reddening 
coefficient were also created for the two maps taking the mean 
error (-0.1) as binning. Figure[3]shows the spatial distribution of 
the two derived c(H/3) maps and their corresponding histograms. 

The structure of the reddening map of the central region is 
mostly uniform with values ranging from 1.3 to 2.5 and a mean 
value of -1.85+0.10; the histogram reveals that the most proba- 
ble value for c(H/3) in this zone is 1.90. Single pixels with very 
high or very low values have a large error in the coefficient esti- 
mations. Big holes in the map correspond to the masked pixels. 

The derived reddening coefficient map of the edge pointing 
has a less homogeneous structure, with a c(H/3) mean value of 
-2.1 1+0.08 over the 1.7-2.8 range. It is interesting to notice that 
all the pixels with c(HJ3)>2.5 are placed in the NW area, where 
the discontinuity of the Ha image is observed (Fig. [2]). In this 
region, pixels with c(H/3)>2.5 were inspected individually; after 
checking the S/N of the Balmer lines and the c(H/3) errors we 
decided not to mask them and to pay special attention to the rest 
of properties derived there. We study this region in detail in the 
ID analysis (Sect. [4|. The statistical frequency distribution of 
the reddening coefficient for the edge pointing gives 2.0 as the 
most probable value. If we exclude values higher than 2.5, the 
distribution can be fitted by a Gaussian function. 

To compare our results with the literature, we estimated 
the e xtinction as A v = 2.145 x c(HJ3), using the |Cardelli et al.| 
(1989) extinction law with R Y = 3.1 and the colour excess as 
E(B - V) = 0.692 x c(HJ3). The mean values derived were 
A v =3.9 and E(B-V)=1.3 for the central pointing, and A v =4.5 
and E(B-V)=1.5 for the edge. The reddening coefficien ts derived 
from our data are higher than those estimated by Esteb an et 



The electron density (n e ) maps were produced from the 
[Sn]/U67 17/6731 ratios using the IRAF packag e TEMDEN 
based on a five-level st atistical equilibrium model ( |De Robertis 
|et al. 1987 ; Shaw & Dufour 1995). Using these maps, we created 
statistical frequency distribution of the electron density with a 
binning of 100 cm -3 (low density limit). They are shown with 
the derived n e maps in Fig. [4] In general terms, the values of 
n e presented in these maps are in good agreement with values 
reported in the literature. 

The histogram for the central pointing shows elements dis- 
tributed in a wide range of density from -200 to -3000 cm -3 
with 1000 cm -3 as the most probable n e . The mean value of the 
distribution is 1008 cm -3 . Some isolated pixels appear very in- 
tense in the image with densities as high as 3000 cm -3 , but with 
large errors. The density distribution follows the low-ionization 
emission elements supporting the idea of a bipolar structure. It is 
interesting to notice that the knots with the higher surface bright- 
ness correspond to the denser zones. 

The histogram for the edge pointing shows a distribution 
close to a Gaussian centred on 500 cm -3 . The map ranges from 
100 cm -3 to 1000 cm -3 with a mean value of 507 cm -3 . As 
happens in the other region, some pixels (7) show higher den- 
sities (up to 1000 cm -3 ), and we removed them from our estima- 
tions. The majority of the pixels with high reddening coefficient 
(c(HJ3)>2.5) were rejected by the S/N mask for the sulphur line, 
but the unmasked pixels present a mean density of 613 cm -3 . 

The morphological analyses showed that the bright knots are 
aligned in a preferred axis along the NE-SW direction with a 
bipolar structure (see Fig. [2]); to check that the electron density 
is related with the bipolarity, we performed a cut in the density 
maps along this direction (see cuts in Fig. [4]), the density pro- 
files obtained are presented in Fig. [5] For them, we performed 
four fits using the least- squares method: the first from the star 
towards the SW, the second from the star towards the NE includ- 
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EDGE [Sii]/Ho 

0.15 

0.12 



CENTRE 



Distonce from the stor (pc) 



Fig. 5. Radial variation in electronic density (in cm -3 ) with dis- 
tance (in pc) along the direction of bipolarity (from NE to SW). 
We consider negative radius from the star towards the NE and 
positive from the star towards the SW. Lines indicate least- 
squares fits: solid lines correspond to pixels differentiating the 
two pointings, and the dashed line represents the fit along the 
direction star-NE including pixels from both pointings. 



ing pixels from the two pointings (centre and edge), and the last 
two fits from the star towards the NE, but differentiating the two 
pointings (see Fig. [5]). It can be seen that the density decreases 
when we move away from the WR star. In addition, the fits show 
a symmetric gradient in the central points with a tendency to 
flatten out towards the ends. 



3.2. Emission line relations 

Figure [6] shows maps of the emission line ratios for the two 
pointings. Their mean values are summarized in Table [3] All 
the intensities presented in the table and figure are reddening- 
corrected. 

In both regions the [Sn]/Li67 17,673 l/Ha map presents an 
inhomogeneus and patchy structure. [Sn] lines are fainter than 
Ha in all the spaxels, with a maximum logarithmic ratio of - 
1.3 to the north of the edge pointing. Some isolated pixels show 
higher values, but they are over the limits of the region masked, 
so unreliable. 

The distribution of the [Nii]^6584/[Sii]^67 17,6731 map 
presents a structure opposite to [Sii]/Ha in both regions. The 
[Nil] emission is stronger than [Sn], reaching log([Nn]/[Sn])=1.4 
in areas close to the ISM. 

Studying the [Nn]/Li6584/Ha maps led to more interesting 
results. The central pointing shows positive values, except in 
some regions in the direction of the bipolarity, where [Nn] and 
Ha fluxes are equal. In the edge pointing, regions with different 
ratios are clearly separated. In most of the pixels, [Nn]>Ha with 
an increasing ratio towards the side. To the north, an area can be 
seen where Ha >[Nn]; this region possesses the higher derived 
c(Hj3), and it was masked in the sulphur maps because of its low 
S/N (<5). The NW area has the highest ratio, possibly produced 
by the contamination of a nearby field star. 

To understand the differences of the [Nn]/L16584/Haf ratio in 
the edge pointing, we complement the study by generating sta- 
tistical frequency distributions of the ratio map and plotting all 
the spaxels from emission line maps in the diagram [Nn]/16584 
vs. Ha (Fig.[7]).The diagram [Nn]/16584 vs. Ha showed in Fig. [7] 
a, presents double behaviour. We considered two lines with unity 




Fig. 6. Derived maps of the emission line ratios of 
the two pointings: edge (left) and centre (right). Top: 
[Sn]^67 17,673 l/Ha. Middle: [Nii]^6584/[Sii]^67 17,6731. 
Bottom: ([Nn]/L16584/Ha. Orientations and sizes as in Fig. [2] 

Table 3. Mean values of the emission line ratio maps 



Line ratios 



Edge Centre 



log([Sn]^67 17,673 l/Ha) -1.03 -1.01 

log([Nn]i6584/[Sn]^6717,6731) 1.15 1.07 
log([Nn]^6584/Ha) 0.07 0.06 



slope as upper and lower limits. Points above the upper line 
are pixels where the [Nn] emission is stronger than Ha, while 
below the lower line they have the opposite behaviour. Points 
between these two lines are pixels with log(Hof)=log([Nn]) ± 
0.05. Then, we located the points of the diagram in the FoV of 
PPAK, taking these limits into account, to identify their spatial 
locations (see Fig. [7] b); they appear grouped. The statistical 
frequency distribution of the [Nn]/Hof map shows a bimodal 
distribution as we can see in Fig(7] c. When we identified the 
spaxels of the three regions defined above, we found that the 
left peak (centred in — 0.3) includes all the points behind the 
lower limit, and the right peak (centred in ~0.1) includes the 
points of the two other zones. We can conclude that at least two 
spatial regions exist in this pointing: one with [Nn]>Hof to the 
SW and another one to the north with [Nn]<Hof. All the pixels 
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Fig. 7. Relations between [Nn] and Ha for the edge pointing. 
Colours help us locate points spatially: red corresponds to points 
with log(Har)=log([Nn]) ± 0.05, blue to points with log(Ha)> 
log([Nn]), and green to points where log(Ha)< log([Nn]). From 
the top to the bottom: (a) log([Nn]/l6584) vs. log(Har). All the 
spaxels of the intensity maps (in units of log(erg cm -2 s" 1 ) are 
represented in the diagram with crosses. Black lines with unitary 
slope represent the limits, (b) PPAK FoV of the edge pointing 
with the zones defined in plot a. (c) Statistical frequency distri- 
butions of the log([Nn]/Hor) map. Black solid line represents the 
distributions of all the spaxels, and coloured dashed lines repre- 
sent the regions defined above. See text for details. 



with c(H/3)>2.5 are included in the second region, along with 
the spectra with very low S/N of the sulphur lines. 

For the central pointing, the same analysis shows that [Nn] 
follows Ha for all the points in a one-to-one relation line of 
unity slope. Relations between the other lines were also stud- 
ied by means of two diagrams ([Nn]/i6584 vs. [Sn]/U67 17,6731 
and [Sn]/U67 17,6731 vs. Ha), showing strong correlations in 
both the pointings. The statistical frequency distribution of all 
the emission line ratios showed single peaks with distributions 
close to Gaussian functions, except [Nii]/Li6584/Hof on the edge. 




Fig. 8. Relative radial velocity field derived for Ha in units of 
km s" 1 . Zero is the global mean velocity (139 km s" 1 ), see 
text for details. The two pointings are mosaicked. The red cross 
marks the position of the central star. North is up and east left 
(see Fig. [T]). 



3.3. The radial velocity field 

Limitations on the instrument resolution prevented us from car- 
rying out an exhaustive analysis of the kinematics of Ml -67. 
Nevertheless, the resolution was sufficient for studying the dis- 
tribution of the radial velocity field and relating it to the mor- 
phology and ionization structure. 

Using the central wavelength of the Gaussian fit performed 
in the cubes, we created radial velocity maps for the two ob- 
served regions. Two corrections were carried out over the mea- 
sured radial velocities. First, we estimated the error in the wave- 
length calibration by comparing the wavelength of a sky emis- 
sion line with its theoretical value, we obtained a difference of 
-0.303 A (—16 km s _1 for [Ol]A55UA% and this zero point was 
added to the measured velocities. Then, we translated maps into 
the local standard of rest (LSR) and corrected for the Earth's mo- 
tions, taking coordinates and universal time of the observations 
into account. 

With the corrected radial velocity fields of Ha, we scaled 
the measured velocities using the overlapping region of the 
two pointings to avoid deviations. Then, we calculated the 
total mean velocity, obtaining a value of 139 km s" 1 , and we 
established it as the heliocentric velocity of the nebula. This 
velocity is in very good agreement with the 137 km s" 1 obtained 
by |Sirianni et al.| ( [l998| ). We present the relative radial velocity 
field of Ha for the two regions mosaicked in Fig. [8] 



Previous kinematic studies (Solf & Carsenty 1982; Sirianni 
et al. 1998) have found two components (one redshifted and an- 
other blueshifted) supporting the idea of a shell in expansion. 
With the low resolution of our data we cannot resolve both com- 
ponents, and the velocity field shown is dominated by the ra- 
dial velocity of the brightest knots, a kind of intensity- weighted 
radial velocity distribution. Despite the low resolution, a study 
of the overall structure of both regions can be carried out. The 
gas of the nebula seems to move faster near to the WR star, de- 
creasing its relative velocity when moving away from the centre. 
The velocity field changes its tendency (increasing) in the "pe- 
culiar" zone towards the north of the edge pointing, where other 
properties were also found to differ from the rest of the nebula. 
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Fig. 9. Statistical frequency distributions of the radial velocity 
of Ha relative to the heliocentric velocity, with a binning of 
5 km s" 1 . The central region is on the top with all the spaxe- 
les represented. On the bottom the edge pointing: black solid 
line represents all the pixels, the short-dashed red line is the re- 
gion where [Nn]>Hof, and the long-dashed blue line represents 
the regions with [Nn]<Ha. 



Figure [9] shows the statistical frequency distributions of the 
radial velocity maps with a binning of 5 km s" 1 . To consolidate 
the differences found in the diagram [Nn]/16584 vs. Ha for the 
edge pointing (Fig. [7|, we represented pixels from the two re- 
gions separately The region where [Nn]>Hof presents a Gaussian 
distribution that covers a wide range in velocity, suggesting that 
some regions are moving away from us and others towards us. 
The distribution of the regions where Ha emission is higher than 
[Nil] is narrower and it is centred near to the zero velocity. 



4. Properties of the integrated spectra 

We created ID spectra by combining fibres to describe the inte- 
grated properties of several interesting zones. Eight integrated 
spectra were generated over the two pointings (see Fig.fTO). 



The regions selected over the central pointing are: 
-Region 1 (Rl): examining the emission of the low ionization 
elements showed in Fig.[2] three bright knots appear to the south 
of the nebula. Eight fibres over these knots were selected and 
combined to create a single spectrum. The offset from central 
star is Aa -4.05", A£ -13.5". 

-Region 2 (R2): we combined three spaxels to the north of 
the star coinciding with another isolated knot. The offset from 
central star is Aa -1.35", AS -14.85". 




Fig. 10. Ha images of the two areas of Ml -67 observed with 
PPAK. Boxes represent the eight regions where the integrated 
spectra were generated. For the offsets of each region from the 
central star (green cross), see the text. Orientations and sizes are 
as in Fig. [2] Edge on the left and centre on the right. 



-Region 3 (R3): we chose those fibres placed to the east of 
the star in a zone where an extended emission is seen in Ha, 
paying attention to not include light from any star. The offset 
from central star is Aa -12.15", AS -4.05". 

-Region 4 (R4): we were interested in analysing a large 
region in the NW masked in the 2D analysis where all the emis- 
sion line maps showed S/N lower than 5. The fourth integrated 
spectrum was created there to check that there is emission in this 
area. The offset from central star is Aa -14.85", AS -14.85". 

The regions selected over the edge pointing are: 

-Region 5 (R5): nine fibres were selected on the south of the 
edge pointing, close to the discontinuity. This spectrum belongs 
to the region showed in Fig.[7]b where [Nn] is stronger than Ha. 
The offset from central star is Aa -31.05", AS -16.2". 

-Region 6 (R6): we combined several spaxels at the SW 
limit of the FoV to check that [Nii]~Hof in this area, as we 
found in Sect. |3.2| The offset from central star is Aa -14.85", 
AS -20.25". 

-Region 7 (R7): in a faint region to the north of the edge 
pointing where interesting properties were obtained in the 2D 
analysis: some pixels show c(HJ3)>2.5, the S/N of the sulphur 
lines is very low so they were masked in several maps, the 
[Nn]/Hof ratio has its minimum values, and the kinematic study 
revealed that here the radial velocity increases opposite to the 
general trend. Seven spaxels were combined in this region to 
analyse the properties in detail. The offset from central star is 
Aa-27",A£-40.5". 

-Region 8 (R8): six fibres were selected on the left of the 
discontinuity to checked whether this region has nebular emis- 
sion. The offset from central star is Aa -47.25", AS -52.65". 

In addition, another three integrated spectra were extracted 
to perform several tests. From the central spaxel of the central 
pointing FoV, we obtained the spectrum of WR124 (Region 
WR). The other two were extracted at -15" to the NE of the star 



(common region i n both pointings) to study the zone where Chu 
|&Treffers|fl98lt found emission in [Om]^5007A (Regions SI 
and S2 in the central and edge pointings, respectively). Figure 
TT] shows six representative ID spectra from the 11 created. 



Fluxes of the main emission lines were measured by fitting 
Gaussian functions using SPLOT within IRAK All the mea- 
sured fluxes are in units of erg cm -2 s" 1 per fibre (area of fibre 
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Fig. 11. Examples of integrated spectra. From left to right and top to bottom: (a) Whole spectrum of Region 3. (b) Spectrum of 
Region 5 in the range of Ha. (c) Region 7, in same range as b, where the absence of sulphur lines can be seen, (d) Spectrum of 
Region 5 centred on the [Nn]/15755A emission line, used to calculate the electron temperature, (e) Zoom over the spectrum in Region 
S2 without any emission in the [Om]/15007A emission line.(f) Whole spectrum of WR124 obtained from the central spaxel. 



-5.7 arcsec 2 ). Statistical errors were estimated using the formula 
presented in Perez-Montero & Diaz (2003]): 

cr l = o~ c N l/2 [l + EW/(NA)] l/2 (1) 

where <xi represents the error on the observed line flux, N is 
the number of pixels used to measure the line, EW the line 
equivalent width, <x c the standard deviation in the continuum 
close to the line of interest, and A represents the dispersion in 
A/pix. 

We derived the reddening coefficient, c(HJ3), from the 
Ha/H/3 and Hy/H/3 line ratios using the procedure described 
in Sect. |3.1| In Region 7, Hy was measured with low S/N, 
and only the other two B aimer lines were used to estimate 
c(Hj3). The reddening coefficients agree with values obtained 
in the 2D study. Table [7] lists the reddening-corrected fluxes of 
the emission lines measured in every zone labelled with their 
standard identification. The third column rep orts the adopted 
reddening curve using extinction law by Card elli et al. ( 1989| ) 
with Ry = 3.1. Errors in the emission line intensities were 
derived by propagating the observational errors in the fluxes and 
the reddening constant uncertainties. The estimated fluxes and 
errors were normalized to F(Hj3) = 100. The values obtained 
for c(Hfi) are also presented in the last row of Table [7] 

Five integrated spectra deserve special attention. First, 
the R4 created in the NW dark area of the central point- 
ing only showed three emission lines: Ha, [Nn]/16548A, and 
[Nn]/i6584A. We deduce that, in areas out of the bipolar struc- 
ture a faint, but not negligible, emission exists coming from the 



nebular gas rather than the ISM. We estimated the Hfi flux by 
means of the reddening coefficient: assuming c(H/3)=1.87 (the 
mean value of the other integrated spectra of this pointing), we 
performed the inverse process of the extinction correction and 
obtained F(#/?)=3.21 10" 16 erg cm" 2 s _1 . 

The spectrum of R8 does not show any emission, thus physi- 
cal and chemical properties could not be estimated here. We did 
not include this region in tables. We extracted a spectrum of the 
WR star (Region WR) shown in FigfTTf, and not perform a de- 
tailed analysis here. Finally, the study performed over Regions 
SI and S2 revealed typical nebular spectra (very similar to R4), 
but we did not find any emission of the [Om]/15007A line as seen 
inFig.[TTfc. 

4.1. Physical properties and chemical abundances 

Electron density (n e ) was calculated from the [Sn]/L167 17/6731 
line ratio using the IRAF package TEMDEN. The derived 
density ranges from -1500 cm -3 near the star, to -650 cm" 3 
towards the edge. These val ues are consis tent with our 2D maps 
and with previous studies (Est eban et al. 1991 Sirianni et al. 
[T998] ). 



^N2 



Electron temperature, T e , can be derived using the line ratio 
7([Nn]^6548) + 7([Nn]^6584) 



7([Nn]/l5755) 



(2) 



The [Nn]/i5755A auroral line that appears close to the "sky" 
line Hgi 5770A, was detected in two zones (R5 and R6). We 
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measured this line again in the spectra before sky subtraction 
and conclude that the flux of [Nn]/i5755A line in R6 is contam- 
inated by the Hgi emission line, thus not reliable. We obtained a 
direct estimate of T e ([Nn]) from R N2 only for R5. 



( |Perez-Montero & Conti ni 2009 ) to estimate N/O from the ni- 
trogen and sulphur emission lines: 



log WO) = 1.26 x N2S2 - 0.86 



where 



To reinforce the validity of the chemical abundances esti- 
mations and to provide ionization correction factors (ICFs) for 
those species whose ionizations stages were not all observed in N2S 2 = log 
the optical spectrum, we performed photoionization models of 
R5. To do so, we used the code CLOUDY v. 10 (Ferland et al. 
1998), assum ing a central i onizing source from a WR star at- 
mosphere ( [Smith et aL||20021 ) with Z = 0.008 and an effective 
temperature of the star of 45 000 K which are, respectively, the 
closest values to the measured total metallicity of the gas and the 
estimated temperature of WR124 ( [Hamann et al.|20 06). 

We considered a spherical geometry putting the gas at a 
distance of 1 pc from the star and assumed a constant density 



7([Nn]26584) 
7([Sn]^6717,6731) 



(3) 



(4) 



Before, we estimated the N/O in R5 with the N2S2 parame- 
ter and checked that the result was in good agreement with the 
value obtained with the direct method. In Table [8] we present all 
the ionic and total abundances, with their corresponding errors, 
deri ved for the integrated spectra. We discuss the results in Sect. 
EH 



of 700 cm -3 , a value similar to the one derived from [Sn] 
emission lines. The model that fits the emission line intensities 
of [On], [Om], Hei, [Nn], and [Sn] better in R5 was obtained by 
varying the ionization parameter (U) and the relative chemical 
abundances of He, O, N, and S. The emission lines from this 
model are listed in Table [7] while the derived physical properties 
and the ionic and total chemical abundances are listed in Table 
U The ICFs obtained were ICF(N + )=1.21 and ICF(S + )=1.58. 
Regarding the resulting geometry, the final radius is 1.22 pc, 
which is of the same order of magnitude as the apparent size of 
the nebula in the images. 

To estimate chemical abundances, electron density and elec- 
tron temperature are required. We used T e ([Nn]) as temperature 
representative of the low ionization ions, S + , N + , and + , and 
T e ([Om]) for deriving the 2+ and He + abundances. In those 
zones where the electron temperature was not calculated, we 
adopted the value of R5. In previous studies, T e ([Nn]) ranges 
from 5900 K ( [Sirianni et aL][T998] ) to 8000 K QChu & Treffers] 
1981); maybe, the supposition of T e =8200 K leads to our 
abundances being underestimated. Since the photionization 
model predicts T e ([Nn])~8550 K and T g ([Om])~8330 K in R5, 
we considered T e ([Nii])~T e ([Oin]) in the estimations. To infer 
abundances in R4 and R7, where sulphur lines were not mea- 
sured, we adopted the electron density of R5, n e =631 cm -3 . We 
checked that variations in density do not affect this estimation. 

Ionic abundances were derived from the forbidden-to- 
hydrogen emission line ratios using the functional forms given 
by Hagele et al.| ( |2008) , which are based on the IRAF package 
IONIC. We used equations from [Olive & Skillman| p004| ) to 
obtain the singly ionized helium abundance. To determine the 
total abundance of O/H we added the two ionic abundances 
(O/H - + /H + + 2+ /H + ). The total N/H and S/H abundances 
were inferred thanks to the ICFs obtained in the photoionization 
model, X/H~(X + /H) x ICF(X + ). In the case of helium abun- 
dances we used the relati on between X(Q 2+ ) =0 2+ /(0 2+ +0 + ) 
and ICF(He + +He ++ ) from jlzotov "eTal] pOOT) Fig. 7) and we 
deduced that ICF(y + )»l. Since our helium measurements 
are uncertain, we do not venture to estimate the total helium 
abundances. 

In R5 all the useful emission lines were measured and the 
abundances determined as explained above. In the rest of the re- 
gions, we did not measure all the necessary lines to calculate 
abundances, and we resorted to the empirical parameter N2S2 



5. Infrared study 

To enhance the morphological and chemical analysis, a study 
in the mid-infra red was performed. We obtained IRS (Infrared 
Spectrograph, |Houck et al. 2004 ) data in mapping mode and 
the MIPS (Multiband Imaging Photometer, |Rieke et aL]|2004| ) 
24 /mi image from the Spitzer Heritage Archive (SHAQM1-67 



van 



has already been studie d in the past in the infrared range by 
der Hucht et al. (1985). They presented the energy distribution 
of the central star WR 124 and flux densities, finding thermal 
emission of dust at T c -100 K. 



Figure 12 shows the MIPS 24/mi image of Ml -67. This 
image has already been presented by |Gvaramadze et aL] ( |20 10 ). 
In a nebula, the origin of the 24/mi emission can be mainly due 
to two factors: presence of the [Oiv]25.90/mi line from highly 
ionized gas or warm dust. Since Ml -67 presents a low degree 
of ionization, we deduce that the observed emission shown in 
Fig. [12] traces the warm dust distribution of the nebula. The 



emission has an elliptical shape along the NE-SW direction 
in very good agreement with the bipolar axis observed in Fig. 
[2j thus suggesting that the structure is composed of a mixture 
of ionized gas and warm dust. Furthermore, an external and 
spherical structure can be seen extending around the ellipsoidal 
shell. This faint bubble is not seen in our optical images. 

For the low-resolution short-low (SL) and long-low (LL) 
modules spectroscopic observations, basic calibrated data 
(BCD, pipeline version 18.18 ) were processed a nd analysed 
with the CUBISM software ( [Smith et al] 12007ft . Data were 
background- subtracted using averaged off-source observations 
and flux-calibrated. Bad pixels were removed with the automatic 
correction routine within CUBISM and a datacube assembled 
for each module. CUBISM allows extracting of ID spectra over 
polygonal apertures: given the different spatial coverage of the 
SL and LL module, we chose two apertures (with an area of 
~60 arcsec 2 ) on the outskirts of the nebula, observed by both 
modules. The spectra from the different modules were stitched 
together, ignoring the noisy region at the red end of each order. 
We called them Regions A and B (see Fig. [12]). In Fig. [13] we 
present the spectrum obtained in Region B. 

We measured the most important lines by fitting Gaussian 
functions with IRAF. Errors were calculated as we explained 
before (Sect. [4]). Fluxes and their corresponding errors are 

2 Website: sha.ipac.caltech.edu/applications/Spitzer/SHA 
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Fig. 12. MIPS 24yum image of Ml -67. North is up and east left. 
Boxes indicate the two regions where IR spectra were obtained. 
Contourns represent the Ha emission derived from Fig.[T] 



(S ++ /H + ) fi . We assumed that S 3+ /H + is negligible. Results are 
presented in Tableland discussed in Sect. 6.1 



Table 4. Lines measured over the two spectra studied in the in- 
frared range. Integrated fluxes are in units of 10" 5 erg cm -2 s" 1 . 





A (jim) 


F(A) 


Line 


Region A 


Region B 


[Siv] 


10.51 




4.1 +0.5 


H(7-6) 


12.37 


4.9 + 0.5 


5.1 + 1.0 


[Nen] 


12.81 


121.7 + 3.2 


105.3 + 3.5 


[Nem] 


15.56 


5.0 + 0.4 


1.1 +0.3 


[Sin] 


18.71 


133.9 + 5.5 


99.2+ 1.9 


[Sin] 


33.48 


156.6 + 5.3 


135.2 + 4.7 



Table 5. Ionic and total chemical abundances estimated in 
Regions A and B with infrared spectroscopy. 

Region A Region B 



[Nell]X12.8 



,WW/' 



[SIII]X18.71 



[SIII]X33.48 




Wavelength (jiim) 



Fig. 13. Infrared spectrum obtained in Region B. The most rele- 
vant lines are indicated. 



presented for the two regions in Table [4] 

Assuming the electron temperature of Region 5 (T e =8200 K) 
and an electron density n e =600 cm -3 , we inferred the chemical 
abundances. To obtain the fluxes relatives to Hj3 we used the 
theore tical ratio of H(7-6)/H/3=0.0109 from |Storey& Hummer 
( |1995| ). The ionic abundances, Ne + /H + , Ne 2 7H + , S z+ /H + , and 
S J+ /H + were inferred by using the IRAF package IONIC. We es- 
timated the total neon abundance by adding the two ionic abun- 
dances, Ne/H~Ne + /H + +Ne 2+ /H + . For deriving the total S/H 
abundance we need to add the S + /H + from the optical spectra. 
To do so we compared the regions from which the IR and the 
optical ID spectra were taken. Noticing the proximity between 
Region A and R3, we approximated the total sulphur abundance 
in the spectrum A as S/H~(S + /H + )tf 3 + (S ++ /H + ) A . The ID spec- 
trum nearest to B is R4, but in R4 we did not measure sulphur 
lines. Since the S + /H + is similar in all the integrated spectra, we 
considered the mean value, (S + /H. + ) mean =6.ll , so that the total 
sulphur abundance in B can be written as S/H~(S + /H + ) mean + 



12+log(Ne + /H + ) 

12+log(Ne 2+ /H + ) 

12+log(S 2+ (18.71//m)/H + ) 

12+log(S 2+ (33.48//m)/H + ) 

12+log(S 3+ /H + ) 

12+log(Ne/H) 
12+log(S/H)t 



7.56 + 0.04 
5.85 + 0.05 
6.59 + 0.09 
6.64 + 0.16 



7.57 + 0.04 
6.72 + 0.07 



7.47 + 0.08 
5.18 + 0.16 
6.44 + 0.17 
6.56 + 0.21 
4.41+0.19 

7.48 + 0.08 
6.63 + 0.11 



t Assuming S + /H + from the optical spectroscopy. 



6. Dicussion 

We included Ml -67 in our IFS observational programme to pro- 
vide answers to some questions that still surround this object: 
degree of gas homogeneity (both kinematic and chemical), stel- 
lar evolutionary phase origin of the gas, interaction with the ISM, 
influence of the star spectral type at WR stage, etcetera. To do 
this, we put together our results for the optical (ID + 2D) and 
infrared analysis, and we complemented them with theoretical 
models of stellar evolution and previous kinematic studies of this 
nebula. 

6.1. Chemical content of Ml -67 

The chemical abundances derived from the ID optical and 
infrared studies, presented in Tables [5] and [8] give us relevant 
information on the chemical content across the nebula. To 
compare the derived abundances with the expected ISM values 
at the location of the nebula, we use the solar values from 
Aspl und et al.| ([2009 ) as our primary reference. For the sake 
of consistency as reference for Ml -67, we consider here gas 
abundances derived following the same methodology, i.e. Hn 
region collisional emission lines. We adopted the chemical 
abundance s of the prototypical Hn region M 42 as a reference 
(Sanchez et al.||2007[|Tsamis et al.||201l| with t 2 =0, and refer- 



ences therein). Then, we corrected the t -0 abundances from 
the effect of the radial abundance gradient of the Milky Way 
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( Rudolph et al. 2006| ) to the galactocentric radius of Ml -6^] We 
considered the constant ratio log (Ne/O)=-0.73±0.08 since they 
are products of the same nucleosynthesis. After these correc- 
tions the expected ISM abundances to be compared with Ml -67 
are 12+log(0/H) -8.42+0.03, 12+log(N/H) -7.54+0.09, 
12+log(S/H) -6.99+0.12, and 12+log(Ne/H) -7.69+0.09. 

First of all, it can be observed that our derived oxygen abun- 
dances in R5 and R6 (12+log(0/H) -7.73+0.06, 7.67+0.07, 
respectively) are substantially lower than the expected value by 
factors -10 with respect to the solar reference, and -7 with re- 
spect to the ISM. This result implies that in the Ml -67 nebula 
oxygen is strongly under-abundant. Comparing the derived N/H 
abundance with the expected ISM value, we find that nitrogen is 
strongly enriched in Ml -67 (factor > 6). 

Overall, this chemical composition can be seen in all the 
nebular regions observed; the N/O ratio appears extremely high 
due to the effect of both nitrogen enhancement and oxygen defi- 
ciency. This fact can be understood when assuming we are see- 
ing regions composed of material processed in the CNO cycle. 
This result for N/O abundance is consistent with previous ID 
studies (Esteba net al.|199l] ), but here it has been extended across 
the whole (2D) nebular geometry and physical conditions. The 
only region where the N/H abundance is close to the ISM ex- 
pected value is R7 (the region with different properties in the 2D 
analysis, see Sect. [3]). 

We did not estimate the total helium abundances since 
our helium lines are very faint and the measures uncertain. 
Nonetheless, given the low limit of the value of Hei (<0.03), the 
absence of Hen and the ICF inferred from Izo tov et aT] ( |2007| ) 
(ICF(y + )»l), we deduced that in Ml-67 the largest part of 
helium is unseen and in neutral form. 

The analysis of the chemical abundances obtained here is 
reinforced by the information derived from the infrared study. 
The infrared spectrum allowed us to derive the sulphur and neon 
abundances for the main ionic species, Ne + , Ne ++ , S ++ , and S +3 . 
The total neon abundance derived in Ml-67 is consistent within 
the errors with the expected ISM abundance for the two aper- 
tures (Table [5]). The noble gas neon is not expected to suffer nu- 
cleosynthetic transformation in the stellar interior, and its abun- 
dance should be preserved. 

In the case of sulphur, the derivation of the total abundance 
requires the contribution of the optical S + to be added to the 
ionic fractions derived from the infrared. Once this approxi- 
mation has been assumed, the total abundance of S/H obtained 
is close to, though still slightly lower than, the expected ISM 
value at the galactocentric distance of Ml-67. Thus we cannot 
rule out the possibility that the nebular material could be 
slightly sulphur-poor: either a certain degree of depletion on 
dust or maybe a nucleosynthetic origin (or both) could be at 
work as reported for some planetary nebulae ( [Henry etal.|2012| ). 



Taking the abundance ratios into account, we can obtain 
clear indications of the excitation degree of the nebula. The 
values N + /N -1 and + /0 ++ >1 from the optical and the derived 
ratios of Ne + /Ne ++ and S ++ /S 3+ from the IR study point to the 
very low ionization degree of the gas in Ml-67. The ionization 
parameter obtained from the photoionization model of R5, 

3 Assumed R G ~10 kpc as the the galacto centric distance of th e rep- 
resentative ISM at the location of Ml-67 jEsteban et al.||1992) and 
taking the distance from the Sun to M42 into account (d~0.414kpc, 
Menten et al.|2007) . 



log (U) = -3.84, is fully consistent with this very low excitation 
observed. 

To provide a summary of the chemical abundances obtained 
across the nebula in the optical and infrared ranges, we have 
grouped regions with similar physical and chemical properties 
(whenever possible). In Table [6] we show the results: <1,2,3> 
represents the average of Rl, R2, and R3, <5,6> the average 
of R5 and R6, and <A,B> the average of zones A and B from 
the IR study. In these cases the corresponding parameters were 
estimated as the mean weighted by the error in each zone. The 
two last columns represent the expected ISM values and solar 
abundances from |Asplund et al.] ( |2009| ), respectively. 



6.2. M 1-67 structure 

Although the first observations of Ml-67 showed a nearly 
spherical shape, the high contrast achieved by coronographic 
studies in the inner regions made a bipolar symmetry clearly 
visible ( Nota et aL]| 1995a ). Owing to the field of view of our 
PPAK observations, we cannot detect this bipolarity; however, 
the narrow-band images from the INT and the interpolated maps 
from PPAK (see Figs. [T] and [2]) show that the bright knots are 
aligned along a preferred axis with "holes" in the perpendicular 
direction. The integrated spectrum of R4, confirms that the 
emission in the holes is very faint (i.e. Hfi was not detected). 
Furthermore, the MIPS image from Spitzer (Fig. 12) also 
reveals the bipolar appearance at 24/mi, suggesting that the 
ionized gas is mixed with warm dust. We emphasize that the 
knots are not only regions with high surface brightness but also 
very dense areas where the [Nn]/Ha and [Nii]/[Sii] ratios show 
the maximum values. 



We support the idea of a preferred axis, but either way, is 
the bipolarity the footprint of an ejection from the star? Looking 
at the radial velocity map of Fig. [8] we can see that velocity 
decreases when we move away from the centre (except in the 
far NE where the gas has peculiar properties, see below). This 
agrees with the studies from Chu & Treffers (1981), who predict 
a faster movement near the star, and Siria nni et al.| ( |1998| ) who 
discuss the idea of a bipolar outflow. The spatial distribution of 
the electron density shows similar behaviour: the mean values 
of maps and integrated spectra of the central pointing are higher 
than at the edge (-1500 cm -3 and -650 cm -3 , respectively). 
Also the electron density decreases along the radial cut seen 
in Fig. [5] with a symmetric gradient in two directions (from 
the centre to NE and to SW) and flattening towards the edges. 
Both analyses lead us to think that the preferred axis is not only 
morphological, but is also the footprint of a mechanism that 
could have expelled material in the past and, later, interacted 
with the ISM diluted and decelerated the gas. 

Leaving aside for a moment the discussion of bipolarity, 
there is another striking morphological feature in this object. 
The IR study at 24//m reveals a spherical bubble surrounding the 
bipolar structure. Kinematic studies from Siria nni et al.| (1998) 
show two different motions in the enviroment of WR124: a 
bipolar outflow and an external spherical hollow shell expanding 
into the ISM. In our narrow-band images from the INT, this 
bubble is not detected possibly because the material is diluted 
in the ISM and very weak in the optical range. A simple sketch 
of the proposed structure of Ml-67 is presented in Figure [14J an 
inner region with bipolar or elliptical shape along the direction 
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Table 6. Summary of inferred properties in Ml -67. 



<1,2,3> 



<5,6> 



<A,B> 



ISVP 



Solar^ 



cQ30) 

n e ([Sn]) (cm" 3 ) 
12+log(0/H) 
12+log(S/H) 
12+log(N/H) 

12+log(Ne/H) 

A(log(N/H)) rf 

AQog(Q/H)) d 



1.87 ±0.01 
1581 +49 



1.87 + 0.01 



6.35 + 0.02 

8.13 + 0.01 8.36 + 0.03 



1.90 + 0.02 
677 + 62 
7.70 + 0.03 
6.40 + 0.02 
8.21 + 0.03 



2.15 + 0.04 



7.92 + 0.03 



0.59 + 0.09 0.82 + 0.09 0.67 + 0. 10 0.38 + 0.09 
-0.72 + 0.04 



8.28 + 0.09 c 
6.69 + 0.04 

7.55 + 0.04 

-0.14 + 0.09 c 



8.42+ 0.03 8.69 + 0.05 

6.99+0.12 7.12 + 0.03 

7.54 + 0.09 7.83 + 0.05 

7.69 + 0.09 7.93 + 0.10 



a Expected ISM abundances a tR r , ~10kpc. 
b Solar abundances from Asplund et al. (2009 1. 



c Estimated assuming log (Ne/O)=-0.73+0.08. 

d Variations with respect to the expected ISM abundance. 




Fig. 14. Sketch showing the structure of Ml -67 around the cen- 
tral star WR124: the bipolar axis along the direction NE-SW and 
the spherical bubble. 



NE-SW surrounded by an external spherical bubble. 

The study of the edge pointing region suggests that the gas in 
the NE possesses different properties from the gas of the bipolar 
outflow. In summary, the properties found for this region are a) 
the largest reddening coefficient of the nebula with c(H/3) >2.5; 
b) the only area where we measure [Nn]<Ha and with the small- 
est N/H abundance estimated, close to the solar neighbourhood 
value; c) increase rather than decrease in the relative radial ve- 
locity; d) absence of [Sn]/L167 17/6731 emission lines; e) the min- 
imum HJ3 flux measured (R7) and lack of [Om]/I5007 or helium 
lines. The presence of these properties puzzles us, but here we 
propose a possible scenario to explain the origin of this region. 
The nitrogen abundance of our "peculiar" region points towards 
material not processed in the CNO cycle (e.g. ISM or MS bub- 
ble) while the morphology (Fig. [T]) and kinematics (Fig. 8]) sug- 
gest that it does not belong to th e bipolar ejection. When looking 
at the bow-shock simulations of |van der Sluys & Ear ners (2003) 
and taking the external IR bubble into account, it is possible that 
the high velocity of the runaway WR124 causes a paraboloid- 
like bow shock around the star sweeping up the surrounding 



medium, so that we are seeing the remaining of this bow-shock 
placed in our line of sight. We should bear in mind that the pecu- 
liar region is spatiall y close to the small reve rsed bow- shock like 
structures found by Grosdidi er et al. ( 2001| ) at the NE periphery 
ofMl-67. 



6.3. M1-67: a consequence of the evolution of the central 
star WR124 

The theory of evolution of massive stars can help us explain 
the observed structure. We compare the stellar parameters from 
the central star in Ml -67 (effective temperature and luminosity 
from Haman n~et al.||2006|) wi th the stellar evolution models 
from STARS (|Eggleton||197Tl |Pols et~aL][T995l |Eldridge & 



Tout||2004| ), |Meynet & Maeder| (|2003| ), and the most recent 



models from Ekstror rTet al.| ( [2012] ) to estimate the initial 
mass of the WR star. Regardless of small discrepancies, all 
the models predict initial mass for WR124 of 60-80 MO (J. 
Toala, private communication). The evolutionary scenario 
for a single massive star with 60 Mo <Mi <90 Mo follows 
the sequence 0-Of/WNL< — >LBV-WNL-WCL-SN ( |Meynet| 
et ai7| |2011| ). After spending a normal life as O stars on the 
main sequence (MS), they evolve towards cool er temperatures 
becoming lumi nous blue variables (LBVs) (Humphrey s & 



Davidson 19941) . These stars undergo extremely strong mass 



loss (up tolO- 3 " 4 MO yr l ) through winds and occasionally 
giant eruptions, and thus peel off parts of their stellar en velope 
to form a small LBV nebulae (LBVN) ( |Notaetal.| 1995b] ). LBV 
stars lose their mass so fast that they rapidly evolve away from 
LBV stage to become WR stars. With an initial mass range 
of 60-80 Mo, we can derive that the central star in Ml -67 
has experienced an LBV phase instead of a red or a yellow 
supergiant phase before becoming a WR star. This idea is in 
good agreement with previous studies of the nature of Ml -67 
based on different observational approaches: Ml -67 is very 
likely the imprint of a previous LBV wind instead of a classical 



red super-giant (RSG) wind-blown nebula (Grosdidier et al. 
T998[|van der Sluys & Lamers|2003) . 



The spectral type of the central star (WN8), tells us that it 
is a "young" Wolf-Rayet and that it has most likely entered the 
WR phase recently. Under this hypothesis, we propose that the 
WR winds have not had enough time to substantially interact 
with the previous nebular material and, therefore, the layers and 
observed features originate in stellar material ejected during the 
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MS and/or LBV phases. Considering a representative expansion 
velocity and the linear size of the nebula, we estimate that the 
ejection happened ~5xl0 4 yr ago. This value is slightly higher 
than the LBV phase duration (~1.3xl0 4 yr, | Garcia - Segura et al.| 
1996) thus supporting the hypothesis that the star has recently 
entered the WR phase. 

Taking the physical sizes and morphologies from an 
hydrodynamical simulations of a 60 Mo star as reference 
(Garcia-Segura et al. 1996), it is possible that the external 
bubble of Ml -67 contains material expelled during the MS 
phase, which is very tenuous in th e optical becaus e of the 
dilutio n with the ISM. |Castor et al] £1975} and |Weaver et al] 
( |1977| ) both built models the derive analytical solutions for 
the dynamic evolution of shock bubbles created by interaction 
between the ISM and the stellar wind in the MS phase. 

Several observational reasons have led us to think that the 
bipolar ejection (or axis of preference) is composed by mate- 
rial ejected during the LBV stage. First, the abundances in the 
knots along this axis show enrichment in nitrogen and deficiency 
in oxygen, a behaviour typical of CNO-processed material in 
phases after the MS stagqj It is common in observations of 
LBVN to find very intense [Nn] emission and absence of [Om] 
flNota et al.|[T995bl |Smith et al.|[T998l |Weis & Duschll|2002] ), 
this being also indicative of a low effective temperature and low 
degree of excitation. Furthermore, many of these nebulae show 
clumpy radial structu res (not multiple shell s) and morphologies 
with preferred axes (Clampi n et al.| [T993). The presence of a 
bipolar ejection in Ml -67 enhances the similarity of the nebula 
to oth er LBVN, which almost all di splay some degree of bipo- 
larity ( |Nota et al.|1995b[|Weis|20011 ). In short, Ml -67 shows the 
general properties of LBV nebulae: linear size, total ionized gas, 
velocity field, IR emission, chemical abundances, line intensi- 
ties, and dy namical characteristics; this clearly poin t s to an LBV 
progenitor ( |Nota et al.||1995b| |Lamers et a"L]|2XJ0T| |Weis|[20TT 



among others). 

The idea of Ml -67 being made up of material ejected during 
the LB V stage was suggested in the past by Siria nni et al. 
(1998) based on the total mass of ionized gas , the expansion 
velocity, and the linear size of the nebula. Also Grosdidier et al. 
( |1998| ) explain the clumpy appearance of Ml -67 by assuming 
the interaction of winds in a previous LBV phase. 

Our study depicts Ml -67 as a nebula with two regions: an 
external spherical bubble with material likely produced during 
the MS and an inner nearly elliptical region along the NE-SW 
direction produced due to an ejection in the LBV phase. We are 
observing a WR nebula with LBVN appearance. 



7. Summary and conclusions 

In this work, we have presented the first integral field spec- 
troscopy study of the ring-nebula Ml -67 around the Wolf-Ray et 
star WR124 in the optical range with PPAK. Two regions of the 
nebula were observed and analysed by means of 2D and ID stud- 
ies. We also obtained and analysed IR spectroscopic data and the 



4 We sh ould bear in mind that models that include rotation (Meynet 
& Maeder 2005 ) and recent observations of O stars in the LMC (Rivero 



Gonzalez et al. 



2012) have revealed that some stars of the MS stage can 



MIPS 24yum image of Ml -67 from Spitzer. In the following, we 
present the main results derived from this work. 

1) We obtain maps from the emission lines that allow us to per- 
form a detailed study of the 2D structure of the nebula: 

• Interpolated maps from the main emission lines show a 
clumpy structure with bright knots aligned along a pre- 
ferred axis in the NE-SW direction. The [Oin]/i5007A 
emission is absent over the whole nebula. 

• The spatial distribution of the reddening coefficient 
maps, c(H/3), presents slight variations between the two 
pointings. In the central region c(Hj3) ranges from 1.3 to 
2.5 with a mean value of -1.85, while in the edge point- 
ing the mean is 2.11, ranging from 1.7 to 2.8. 

• Electron density maps, n e , derived from the 
[Sn]/L167 17/6731 ratios, show a non-uniform struc- 
ture. Knots with higher surface brightness in Ha possess 
the highest densities. We also find that density decreases 
with increasing the distance from the star showing a 
symmetric gradient. 

• We analysed the ionization structure by means of line 
ratios maps. In particular, the [Nii]/Ha map of the edge 
pointing field reveals two behaviours, thus defining two 
spatially well delimited regions: one in the NE with 
[Nn]<Hof and the second one with [Nn]>Hof. 

• With radial velocity maps we studied the kinematics of 
the nebula. The derived heliocentric velocity for Ml -67 
is -139 km s" 1 , in agreement with previous results. The 
relative radial velocity seems to decrease as it moves 
away from the central star along the preferred axis. 

2) We derived the physical parameters and chemical abun- 
dances of Ml -67 using the integrated spectra of eight re- 
gions: 

• The electron densities inferred on the central region 
present higher values than on the edge (-1500 cm -3 and 
-650 cm -3 , respectively). This result agrees with the ra- 
dial variations of the 2D study. 

• We derived an electron temperature of -8200 K in R5 
by using our measurement of the [Nn]/15755A emission 
line. 

• The chemical abundances show, in all the studied areas, 
an enrichment in nitrogen and a deficiency in oxygen. 
The nitrogen enhancement in each region is different, 
suggesting an inhomogeneous chemical enrichment. 

3) The 24/mi image reveals an inner bipolar-like structure in 
the NE-SW direction and an outer faint spherical bubble in- 
teracting with the surrounding ISM. From the low-resolution 
mid-IR spectroscopic data, we measured the main emission 
lines and estimate ionic and total chemical abundances, 
verifying the low ionization degree of the gas. 

4) Overall, this study revealed the clumpy structure of Ml- 
67 with knots aligned along a preferred axis and with 
"holes" along the perpendicular direction. The gas along 
this bipolar axis possesses a low ionization degree, and 
it is well mixed with warm dust. The optical analysis 
of these knots revealed chemical abundances typical of 
material processed in the CNO cycle, suggesting that 
the material comes from an evolved stage of the star. 
The radial variations in electron density and velocity indi- 
cate that the gas of the bipolar feature was ejected by the star. 



also show CNO-processed material. 
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5) A region placed to the NE of the nebula shows different 
kinematic, chemical, and morphological properties. We 
propose that this region comprises the remaining of a bow- 
shock caused by the runaway WR124 with ISM material 
mixed up with the MS bubble. 



Based on our observational results and taking theoretical 
models fro m the literature into account (e.g. Garcia-Segura 
et al. 1996| ), we propose a scenario where the central star has 
recently entered the WR phase. This implies that the interaction 
of WR winds with previous surrounding material is not visible 
yet. After comparing our results with stellar evolution models 
and taking the inferred initial mass of the star (60 Mo < 
M; < 80 Mo) into account, we deduced that the central star 
experienced an LBV stage before becoming a WR. The bipolar 
material observed belongs to an ejection during the LBV stage 
since the morphology, kinematics, and chemistry are in good 
agreement with previous studies of LBV nebulae. 
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